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ABSTRACT 

We investigate the propagation of electromagnetic waves in the core of a winding 
number one string defect with isotropic core in nematic liquid crystals. We numerically 
solve wave equations for the TE mode and show the existence of guided modes arising 
due to the variation of the refractive index as a function of the scalar order parameter. 
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1. INTRODUCTION 



Optical fiber waveguides have been a subject of intense theoretical and experimen- 
tal investigations due to their technological importance. Propagation of electromagnetic 
waves in such waveguides has been analysed in detail for different waveguide geometries 
as well as for different optical properties of the materials used. The simplest optical fiber 
essentially consists of an inner cylindrical core of a transparent medium with refractive 
index n\ surrounded by an outer cladding of another transparent material with refractive 
index n 2 such that n\ > n 2 [|]]. There are also graded index fibers where the refractive 
index decreases continuously in radial direction Recently it has been suggested that 
the propagation of electromagnetic waves in the cores of string defects in certain trans- 
parent condensed matter systems may show interesting behavior and that these strings 
may behave as graded optical fiber waveguides @. However, the considerations in || 
were based on geometrical optics notions which are not very appropriate for the realistic 
cases when the thickness of the string defects is of the order of the wavelength of light 
used. 

In this paper we analyze in detail the propagation of electromagnetic waves in the core 
of a winding number one string defects with isotropic cores in nematic liquid crystals. 
Optical fibers with liquid crystalline cores have been discussed in the literature by Lin et 
al. [f|, |5| and we follow the approach used in Q in our investigation. Our work differs 
from that in ||J] in two important aspects. The propagation of light was considered in 
|H for the case when the liquid crystalline core was entirely in the nematic phase (apart 
from the consideration of the effects of the propogating fields on the order parameter 
as discussed in 0), even though the configuration of the director field considered in jj| 
corresponded to strength one line defect. The reason being that for physically interesting 
wavelengths, the cylindrical cavity had to be large making it energetically unfavorable for 
the strength one defect to develop an isotropic core. We, however, consider the case when 
the core of the string is in the isotropic phase. This leads to radially decreasing value of 
the refractive index for the TE mode (as we will see later) resulting in the existence of 
guided modes even without a transparent cladding. This is the second difference between 
our work and that in Q where guided modes were due to the lower value of the refractive 
index in the cladding material. Even though our consideration of the isotropic core for 
the string restricts us to a wavelength regime which is probably not much of practical 
use (in the sense that absorption of light may be strong for such short wavelengths), it 
illustrates a novel possibility where the guided modes arise purely due to the variation 
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of the order parameter field. As the order parameter field varies continuously, we expect 
the refractive index to do the same (as we will explain below) making these strings into 
graded optical fibers. 



2. WAVE PROPAGATION 

We consider a cylindrical waveguide consisting of nematic liquid crystalline core at a 
temperature below the nematic-isotropic transition temperature. The walls of the cylinder 
(which need not be transparent) are assumed to be coated with appropriate material 
so that the director is anchored normal to the wall leading to a strength one string 
configuration |J. We further assume that the radius of the cylinder is small enough so 
that at the center of the cylinder the order parameter vanishes leading to an isotropic 
core for the string. The conditions under which this is energetically favorable have been 
analysed in |4| and we refer the reader to |4[] for further details. The director then lies in 
the plane perpendicular to the axis of the cylinder and points radially from the center. 
We will consider the wave propagation only near the core of the string and hence neglect 
any possible surface effects. 

The radial variation of the order parameter S for a strength one string defect (in the 
single constant approximation) is governed by the Landau-de Gennes free energy density 

F = -S 2 S 3 + — 5 4 + — t + -(-*-) (1) 

2 3 4 2r 2 2 v <9r v ; 

Here we have omitted the surface contribution as we are only interested in the idealized 
case of the wave propagation near the core of the string. Here A = a(T — T ), T is the 
temperature and a, To, B, C and L are parameters of the material. We take the values 
of these parameters [|, as, a = 0.1319 x 10 5 J/m i K, B = 1.836 x 10 5 J/m 3 , C = 
4.050 x 10 5 J/m 3 , L = 2.7 x 10" 11 J/m and T = 307K. With these parameters, the 
string profile can be obtained from the following equation corresponding to the above free 
energy 

+ - 4 - U AS ~ BS2 + CS ") = (2) 
ar z r ar r z L 

We have solved the above equation numerically using a Runge-Kutta algorithm of 
fourth order accuracy. The boundary conditions for the order parameter are that S(r) 
vanishes at r = (corresponding to the isotropic phase) and at large r, S(r) approaches 
its asymptotic value 7] corresponding to the nematic phase. The profiles of S(r) we obtain 
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for various temperatures are in agreement with those in ||. We now try to determine 
the spatial variation of the refractive index for this configuration and whether it can lead 
to the existence of guided modes. At the center of the string the refractive index should 
have the value rii appropriate for the isotropic phase. At large r, in the nematic phase, 
there are two indices of refraction; n» for the extra-ordinary ray and n±_ for the ordinary 
ray. We will consider the propagation of TE mode for which the n± is the appropriate 
refractive index for the nematic phase. As we have mentioned earlier, guided TE modes 
arise here due to the fact that the refractive index decreases from the value n« at the 
center of the string to a smaller value n± away from the center. 

We now assume that the dielectric permittivity e± for the ordinary ray is proportional 
to the order parameter S 0. With the boundary conditions that e±(S = 0) = nje and 
e±(S — 7]) — n 2 j_e , we can derive the spatial dependence of e_i_(r) resulting from the spatial 
dependence of S(r), 

*±(r) = eoK 2 - ^^S(r)) (3) 
where €q is the permittivity of the vacuum. 

To study the propagation of electromagnetic waves in the waveguide we consider the 
Maxwell's equations 

VXE = -f (4) 

VXH = § (5) 

V*D = (6) 

V-H = (7) 

Where, 

D = eE (8) 

B = // H (9) 

here /io is the permeability of the vacuum (we will assume this to be the same for the 
medium as well) and e is the permittivity of the medium. 
We consider the propagating modes in the following form 

E(vi) = E(r)e i(wMz) (10) 
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with similar (z, t) dependence for H as well, uo and (3 are the optical frequency and 
propagation constants respectively. Uncoupled equations for field components can be 
obtained From Eqns.(4)-(9) for the TE mode by taking E = (0, E h 0) and H = (H r , 0, H z ). 
We assume that the director varies uniformly with the azimuthal angle. Since E$ is always 
normal to the director, the permittivity only depends on r and is given by Eqn.(3). 
We get, 

-9/ + " V + [" VoK - - - P 2 - ^}E, = (11) 

H r = — —Es (12) 

UJfiQ 

H. = + (13) 

ufi or r 

Here we have replaced e by e± from Eqn.(3) as appropriate for the TE mode. Using 
eofiouj 2 = (2ir/\) 2 and writing (3 as 2n/\ z , we can rewrite the above equation as 

^ + ;^ + I-4-bs W -^ = o (14) 

B = ^K 2 "<] (16) 

We may mention here that guided TM modes will not be expected to arise for the 
case under consideration. This is because for such modes E r does not vanish and for E r 
the corresponding refractive index ray in the nematic phase, away from the string center, 
is larger than the value rii at the center. 



3. NUMERICAL SOLUTIONS 

We solve Eqn.(2) for various values of temperatures and use the resulting profile of 
S(r) in the expression of e_i_(r) to solve for E<p in Eqn.(14). Numerical solutions are 
again obtained by using a Runge-Kutta algorithm of fourth order accuracy. Boundary 
conditions for E^ are that = at r = and at r — > 00 (for guided mode). We choose a 
fixed value of B and vary A to find the solution of Eqn. (14) which is well behaved at large 
r. We then decrease the value of B to search for the largest value of the wavelength A for 
which a guided mode can be found. Values of A and \ z can be obtained from the values 
of A and B from Eqns. (15)-(16). Profiles of H r and H z are obtained from Eqns(12)-(13). 
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We take the values of the indices of refraction as, rii = 1.6 and n± — 1.5 It is 
important to mention here that these values of refractive indices are typical at optical 
wavelengths while the wavelengths relevent for our case come out to be much shorter. 
Though, we note that numerically we obtain the value of A and B and the refractive 
indices are used only to determine the corresponding values of A and X z from Eqns.(15)- 
(16). We take the above values of refractive indices as an example. Thus, if the value of 
n i ~ n i i s larger for shorter wavelengths then the corresponding wavelength will also be 
larger. [As we mentioned earlier, for short wavelengths the absorption of light may also 
be strong. We will, however, not worry about this. Our main aim here is to illustrate the 
qualitative feature of the existence of guided modes purely due to the spatial variation of 
the order parameter.] 

Fig. (1) shows the plots of and H z for different temperatures and for different 
values of A. H r can be simply obtained from Eqn.(13). Our choice of temperatures was 
guided by the consideration of finding largest values of the core radius of the string which 
would lead to large, possibly practically more interesting, values for the wavelength A. 
The most favorable value of A we could find for our set of parameters is about 2670 A 
corresponding to the plot shown in Fig. la. 

5. CONCLUSIONS 

We have shown the existence of guided modes for a strength one string configuration 
with isotropic core due to the variation of the refractive index as a function of the order 
parameter. Even though the consideration of isotropic core has restricted as to the con- 
sideration of wavelengths of the order of 2600 A, we believe that our results point out 
an interesting possibility that for suitable liquid crystal samples for which the string size 
can be made larger, guided modes can be found purely within the liquid crystalline core 
as compared to the cases where the lower refractive index of the cladding is responsible 
for guided modes. Further, due to continuous dependence of the refractive index on the 
order parameter, the resulting waveguide is similar to a graded index fiber. 
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FIGURE CAPTIONS 



(1) Profile of the order parameter S(r) is shown by the solid curve while and \H Z 
are shown by the dashed and the dotted curves respectively, (a) Plots for AT = 1.44 
and A = 2677A A 2 is found to be = 1790 A. (b) Plots for AT = 1.35, A = 1420l and 
A 2 = 9451. 
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